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We study non-equilibrium current and occupation probabilities of a two-orbital quantum dot. The 
couplings to the leads are allowed to be asymmetric and orbital dependent as it is generically the case 
in transport experiments on molecules and nanowires. Starting from a two-orbital Anderson model, 
we perform a generalised Schrieffer- Wolff transformation to derive an effective Kondo model. This 
generates an orbital potential scattering contribution which is of the same order as the spin exchange 
interaction. In a first perturbative analysis we identify a regime of negative differential conductance 
and a cascade resonance in the presence of an external magnetic field, which both originate from 
the non-equilibrium occupation of the orbitals. We then study the logarithmic enhancement of 
these signatures by means of a renormalisation-group treatment. We find that the orbital potential 
scattering qualitatively changes the renormalisation of the spin exchange couplings and strongly 
affects the differential conductance for asymmetric couplings. 

PACS numbers: 05.60.Gg, 73.63.Kv, 05.10.Cc 



I. INTRODUCTION 

Over the past decade the fast progress in nanofabri- 
cation allowed the observation of single electron tunnel- 
ing as well as Coulomb blockade, cotunneling and Kondo 
signatures in quantum dot devices. Initially these dots 
were gate-controlled isolated islands in two-dimensional 
electron gases like GaAs [l|. Recently a new trend has 
evolved towards studying molecules in gated junctions. 
Starting from small molecules @, Q and bucky balls [4[ 
also carbon nanotubes [fj|, |g] and InAs nanowires were 
studied [7[. In contrast to gate-defined quantum dots 
these molecules and nanowires possess a rich spectrum 
of internal degrees of freedom, which significantly influ- 
ence non-equilibrium transport through these systems. 
In particular, the overlaps of the orbitals of a molecule 
with the states in the leads vary significantly, which re- 
sults in asymmetric and orbital-dependent tunneling am- 
plitudes. 

Parallel to the experimental progress there have been 
numerous theoretical studies on non-equilibrium trans- 
port through nanostructures [8[ . One of the main chal- 
lenges consists in the new energy scale introduced by the 
bias voltage, which is in general of the same order as the 
internal energy scales of the system. The understanding 
of non-equilibrium phenomena in the presence of strong 
correlations represents an outstanding task in condensed 
matter physics which requires the development of pow- 
erful theoretical tools. Recent advances include various 
renormalisation group (RG) methods |fj|-(ll| introduced 
for the study of quantum dots mainly dominated by spin 
fluctuations. 

Motivated by the experimental developments 0t3| we 
will study a quantum dot possessing multiple orbitals 
and asymmetric tunneling amplitudes. A previous per- 



turbative study [121 ] showed that such systems generically 
exhibit interesting non-equilibrium features like the so- 
called cascade effect and the appearance of a negative 
differential conductance (NDC). Our main focus is the 
analysis of the logarithmic enhancements of these fea- 
tures due to the Kondo effect and the interplay of spin 
and orbital degrees of freedom. Even without spin-orbit 
interaction the renormalisation of spin and orbital scat- 
tering processes can influence each other in multi-orbital 
quantum dots. The minimal model to observe these fea- 
tures is provided by a two-orbital quantum dot with finite 
orbital splitting, an external magnetic field, and differing 
tunneling amplitudes (see fig. H]). 

We consider the Coulomb blockade region with a sin- 
gle electron on the dot. In this situation the orbital or 
spin degree of freedom can be changed by elastic or in- 
elastic cotunneling processes. Using a Schrieffer- Wolff 
(SW) transformation [13| we derive a Kondo spin ex- 
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FIG. 1: (Colour online) Sketch of the two-orbital Anderson 
model. The energies of the two orbitals are split by 8. The 
orbitals are coupled to free fermionic leads held at different 
chemical potentials ±V/2. The system is assumed to be in the 
regime of single occupancy and the dot orbitals are subject 
to an external magnetic field B. 



change as well as an orbital potential scattering (PS) in- 
teraction, which have the same strength and therefore 
must be treated on an equal footing. The orbital PS 
term includes two processes, one involving a change of 
the orbital index and one describing transport via the 
empty orbital. 

We first study the orbital occupations and the differ- 
ential conductance to second order in perturbation the- 
ory (PT), which provides a reliable and comprehensive 
picture of the considered non-equilibrium physics. In 
the main part of the paper we perform an RG analysis 
with frequency-dependent couplings [9( and study the ef- 
fect of the orbital PS. This method was recently used to 
study cotunneling transport through carbon nanotubes 
modelled by a two-orbital Anderson dot with two elec- 
trons [6(. While for this situation orbital PS contribu- 
tions are expected to play an important role, the impact 
of an empty orbital remains to be assessed. Furthermore, 
transport through a double quantum dot coupled to four 
different leads was investigated in ref. [14]]. In this work 
NDC with respect to one pair of the leads was observed 
and a detailed RG analysis of the transport properties as 
a function of the orbital splitting was performed. 

In our set-up (see hg. [1]) we study the dependence of 
the cascade resonance and the NDC on both the coupling 
asymmetry and the external magnetic field. By includ- 
ing the logarithmic corrections to the inelastic cotunnel- 
ing processes we further show that the interplay between 
spin and orbital scattering is most pronounced for strong 
asymmetric couplings to the reservoirs. The detected en- 
hancement of the NDC reveals the importance of Kondo 
correlations for the appearance of this effect, whose sig- 
natures can be probed in transport experiments [2j. 



II. MODEL AND SW TRANSFORMATION 

We start from a two-orbital Anderson model sketched 
in fig. [1] The Hamiltonian is given by H = -ffi ca ds + 
Hdot + H lnt . In second quantisation the fermionic leads 
are described by Hieads = J2 a ka( £ aka-fJ. a )c aka c ak a with 
a = L.R, the chemical potentials fi a = ±eV/2 and the 
density of states N(u) = N O(D - |w|), N = 1/(2D ). 
The dot Hamiltonian reads 



n as = d\ s d as . We use a generalised SW transformation 
to derive an effective Kondo model i7x + -ffps in the limit 
8 <C U. We find the conventional Kondo spin exchange 
interaction in its two-orbital form 
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with JP Q = 2t 0a t* ab (l/{e b + U) - l/e„), S ba 

lY,ss> d l s Tss>das>, and Sf3 a = ^J2a*' c l^<T'C a *>- In an 
experimental set-up the overlaps between the states on 
the dot and in the leads will strongly depend on the lead 
and orbital indices, thus leading to a strong asymmetry 
as well as orbital dependence of the hopping amplitudes 
t aa and the Kondo couplings J b " ■ In addition to © the 
SW transformation generates a spin-independent orbital 
PS term 
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where a — 2(1) if a = 1(2). The second term ex Jg" 
describes cotunneling transport across the empty orbital, 
e.g. orbital 2 if orbital 1 is occupied. The first term 
oc Jg" is a spin-conserving process which changes the 
orbital state of the dot; we therefore expect a significant 
influence on the current at voltages V ~ S. We emphasise 
that §5$ is of the same order as the Kondo interaction @ 
and hence cannot be neglected. In particular, we show 
below that the orbital PS term qualitatively changes the 
RG flow of the Kondo couplings. In order to distinguish 
the Kondo and orbital PS contributions we will denote 
the couplings of the latter by P° (= J§£) and P x (= jf") 
in the following. Furthermore, we study the system at the 
particle-hole symmetric point, i.e. s\ — — U/2 — 5/2. In 
this case the conventional single-orbital PS contribution 
(oc l/(e Q + U) + l/e a ) can be neglected ({ < f/ has 
already been assumed in the SW transformation). 

To simplify the notations we introduce a generalised 
vertex V via 



H K + H PS = £ J2 V&iTdldaJ^c^, (4) 
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with the orbital index a — 1,2, the spin index s = ±1, 
and an orbital-independent local Coulomb interaction. 
The energies of the orbitals are £\ s = e\ — ^sg^sB and 
e 2s = £is + 8 with an external magnetic field B. We 
stress that the orbital splitting 6 is assumed to be fi- 
nite. The orbitals are coupled to the leads by H mt = 
H aaa {taaC a<7 d a<7 + h.c.) with spin-independent hopping 
amplitudes t aa and c\ a = £) fe c f akcr . 

In the following we will consider the regime with ex- 



actly one electron on the dot, i.e. J2as 
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Initially both interaction terms © and (O are rotational 
symmetric in spin space, i.e. J z = J 1 - = P° = P x are 
independent of the spin index s on the dot. As we will sec 
below, the RG flow leads to a breaking of this symmetry. 
Furthermore, from now on we use the Einstein summa- 
tion convention and parametrise the couplings as J b " = 
TabrpaJa with r ab = t ab /t L1 and t L1 = t L2 = 1, t m = r, 
tR2 — 1/?'- We use 5 as unit of energy (6 = 1), r as 
asymmetry parameter and set h = ks = e = g/is = 1 • 



III. PT AND RG EQUATIONS 

We first study the system using second order PT. 
The current is computed using a rate-equation approach 
(see e.g. ref. [HI]). Following ref. [l5| we use a pseud- 
ofermion representation [16[ for the Kondo spins; for de- 
tails of the calculation we refer to ref. [17| • We define the 
dot Green's function, G as (T,T') = -i(Tcd as {T)d\ s (T')) , 
where Tc is the contour ordering of the times r, r' on 
the Keldysh contour. For the lesser Green's function we 
assume G^ s (ui) = in as A as (uj), where A as (uj) denotes the 
non-interacting spectral function given by a (5-function at 
e as - The stationary occupation of the state as is deter- 
mined by the rate equation 
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The rates are defined via the retarded self energy as 
J2 bs > Tasks' = r as (e as ) = -2Im[S^ s (e os )]. In second 
order PT (Fermi's golden rule) we find 
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Here w(x) = xub{x), with the finite-temperature Bose 
function ub{x), gives the energy window allowed for tran- 
sitions. We note that the self-consistency equations ([6]) 
have to be solved under the constraint ^2 as n as = 1. 
Results of the perturbative treatment of the occupation 
numbers and the differential conductance are shown in 
figs. 02 and |U 

In order to include Kondo correlations we use a non- 
equilibrium poor man's scaling approach with frequency- 
dependent coupling functions @. The underlying ap- 
proximations are devised for max (V, T, B) >■ Tk provid- 
ing a cutoff for the flow to strong coupling. We start with 
the representation (QJ. During the RG flow the vertex 
V acquires a dependence on the energies of the incom- 
ing and outgoing lead electrons as well as the involved 
dot states. Approximating the latter by their respective 
resonance energies e as and using energy conservation we 
deduce that V will depend on the frequency ui of the in- 
coming electron only. Thus during the RG flow we can 
use the parametrisation ([3]) with w-dependent couplings 
J z , J ± , P° and P x . Explicitly the RG equation for the 
vertex reads 

VLi I. [ui) = 
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where W = Q(D—\w + iT\), D < Dr, is the running cut- 
off and r is the decoherence rate cutting off the RG flow 
9]. Motivated by the expressions entering the computa- 
tion of the current we use the mean value of the rates, 



T = jJ2 as Fas- Generalising the rate ([7]) from PT using 
frequency-dependent coupling functions, we find 
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with the Fermi function f a (uj) = l/{l + e^-^l T ). Here 
we have used that the spectral function A as (uj) is given by 
a Lorentzian of width F centred around ui = s as (replac- 
ing the (5-function) and thus is strongly peaked. The con- 
volution of the Fermi functions l with broadened spectral 
functions leads to an effective temperature cut ting off I 
the flow to strong coupling at low bias voltages [15| . The 
equations © and ([9]) have to be solved self-consistently; 
we find a good convergence within a few iterations. 

Typical results for the flow of the frequency-dependent 
Kondo couplings J are shown in fig. [21 where we focus on 
one representative component of J z ' ± and P°' x . In order 
to assess the importance of the orbital PS contributions, 
a comparison between results with and without orbital 
PS is provided. Lowering the cutoff down to the self- 
consistently determined rate T leads to an overall increase 
of J with respect to the bare initial couplings Jo. Fur- 
thermore, pronounced peaks develop in correspondence 
of resonant transport involving spin- and orbital-flip pro- 
cesses. These peaks are split due to the finite bias by 
±V/2. The singular behaviour is cut off by T; the in- 
duced noise due to the stationary current leads to a finite 
life time of the dot states [9j . This behaviour as well as 
the position of the resonances can be inferred from the 
analytic structure of the flow equations. In the following 
we focus on the coupling J z , the RG equations of the 
other couplings can be found in ref. [17]. Explicitly, the 
RG equation for one representative component of J z in 
the absence of orbital PS reads 
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Characteristic Kondo divergences appear at u = ±(5+B) 
and uj = ±B, which are split by ±V/2 in presence 
of a finite bias voltage. Similarly to the conventional 
Kondo model, J z is renormalised by two spin-flip pro- 
cesses ex (J 1 -) 2 . Schematically, the RG equation has the 
form dJ/(d\nD) = —2 J 2 thus leading to a divergence at 
low energies. 




FIG. 2: (Colour online) Flow of the couplings J z (a), J 
(b) and P°' x (c) as a function of cj for r — 1, B = 0.3, 
T = 0.01, and V = 0.1; the initial values are iVo Jo = 0.01 and 
Do = 1000. Results with (red) and without (blue) orbital 
PS are compared. The evolution of J z under the RG flow 
is indicated by dashed red lines; the final value is obtained 
at D = T w 0.0002. We note that the RG flow results in 
J z 7^ J and that the orbital PS yields new structures in J z 
at w = -5 ± V/2 and oj = 5 ± V/2. The orbital PS is only 
weakly renormalised. 



Qualitatively new effects arise in the presence of orbital 
PS contributions, which lead to the following additional 
terms on the right-hand side of (fTU)l 
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(ii) 

P x always involves a change of the orbital index with en- 
ergy ±8. As can be seen in fig. [2]; the orbital PS couplings 
are only weakly renormalised. For the symmetric case 
the orbital PS couplings P do not flow at lowest order, 
they only evolve due to the finite frequency-dependence 
acquired by the Kondo couplings during the flow. More- 
over, the averages of P x entering the expressions for the 
rates (|9]) and the current (fT2j) will only slightly differ from 
the bare values and thus induce only minor effects. How- 
ever, the flow of P feeds back into the flow of the Kondo 
couplings where it results in qualitatively new features. 



As shown in fig. [2^,, for u> = S ± V/2 the flow of J z (o-i) 
exhibits a crossover from a gradual increase to a dip as T 
decreases. In this regime the flow equation assumes the 
form dJ/(d\nD) = JP/2; the reversed sign with respect 
to the characteristic Kondo divergence gives rise to the 
observed dips. These are not present for B — as they 
merge with the larger peak coming from 5 + B. In this 
case the orbital PS terms only induce an asymmetry as a 
function of u. This can be observed for J- 1 (see fig. [2p), 
where the contribution of the orbital PS only changes the 
weight of already resonant processes. J 1 - is renormalised 
by J^J Z and J P, where both J z and P do not involve 
spin-flip processes and thus no difference appears at fi- 
nite magnetic field. For asymmetric couplings r > 1 we 
find that existing peaks at w = 0, ±i? are enhanced while 
all others are strongly suppressed. 



IV. CURRENT AND DOT OCCUPATION 

The current is obtained by taking into account the gen- 
erated frequency dependence of the couplings. Specifi- 
cally we find 
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with A(lo) a Lorentzian[19j of width 2r centred around 
£bs' — £as- Given the complicated form of the RG 
equations a detailed analytical understanding of the be- 
haviour of the physical observables is not feasible, espe- 
cially as weighted integrals over the coupling functions 
are involved. Thus we have evaluated the dot occupation 
numbers, the decoherence rates and the current numer- 
ically. The renormalisation of the energy levels and the 
magnetic field are neglected in the present analysis, since 
they turn out to contribute only at higher order. 

First we present results for the symmetric situation 
r = 1 and discuss the effects of a finite magnetic field. 
Fig. [5Ji shows the differential conductance obtained from 
PT and RG calculations both with and without orbital 
PS. As long as the couplings are symmetric, the differen- 
tial conductance is symmetric with respect to the voltage 
and thus only V > is shown. We observe that the in- 
clusion of the orbital PS leads to an enhancement of the 
conductance, which is simply due to the increase of trans- 
port channels. This enhancement is voltage independent 
except that its value changes at V = 5, which indicates 
the importance of the coupling of the orbital PS terms 
to the J z processes. Regarding the effect of the RG we 
observe an overall increase of the conductance, which is 
due to the increase of the Kondo couplings under the RG 
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FIG. 3: (Colour online) Differential conductance, G = dl/dV, 
normalised to Go — 2e 2 /h, for symmetric couplings (r = 1); 
a) RG results (full lines) are compared to PT (dashed lines), 
both with (red) and without (blue) orbital PS for parameters 
as in fig. [2] b) RG results including orbital PS for the same pa- 
rameters but different values of B. The cascade step is visible 
at V — S — B provided B < 8/2; Inset: schematic representa- 
tion of the possible transitions between the dot orbitals, the 
cascade transition being marked in red. 



flow (see fig. [5]). Furthermore, there is a clear enhance- 
ment of the inelastic cotunneling signatures at V = B 
and S + B involving spin-flip processes. In contrast, the 
effect of the orbital PS remains unchanged under the RG 
flow. 

The cascade step [12| at V = 5 — B has its physical ori- 
gin in the non-equilibrium dot occupations (see fig. Hk) . 
For V > B the orbital 1J, gets populated from which 
the next higher orbital 2f is split by S — B (see inset of 
fig- Eh)- Therefore, a finite population of 2f already sets 
in at V = S — B, which results in the cascade step in the 
conductance displayed in fig. [3J This is in analogy to the 
Balmer series observed in the traditional spectroscopy of 
hydrogen, where instead of transitions from the ground 
state (Lyman series) the transitions from an excited state 
are probed. As one can observe in fig. [3^ the cascade res- 
onance is only weakly enhanced by RG. The reason might 
be that the RG only starts to play a role at higher or- 
ders since the cascade effect is due to transitions from an 
excited state. Clearly the cascade effect is not present 
at B = while for B = 6/2 it is superimposed by the 
inelastic cotunneling process at V = B (see fig. ^p) . 

In the following we discuss the effect of asymmetric 
couplings and the resulting NDC. Fig. U shows PT and 
RG results for the dot occupations and differential con- 
ductance for asymmetric couplings with r — 5. This 
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FIG. 4: (Colour online) Dot occupation numbers (a) and 
differential conductance (b) for asymmetric couplings with 
r — 5, NqJq — 0.005 and other parameters as in fig. [2] For 
the occupation numbers RG results (full lines) are compared 
to PT (dashed lines) both including orbital PS. 



implies that tm 3> £_R2, i.e. the transport through the 
second orbital is strongly suppressed. At V ~ 6 we ob- 
serve a strong increase of ni-\ at the cost of nif which 
consequently leads to a decrease in the current and thus 
the appearance of NDC. The increase of 7124. similarly 
yields a side peak located at V « 5 + B. A simple expla- 
nation for the appearance of the NDC is provided by PT 
in the limit of r 3> 1, B = and T = 0. Introducing the 
orbital polarisation M\i — ^J2a( Ulcr ~ n 2a) we obtain 
for the differential conductance a.t V ~ 5 



^oc2 + 3M 12 + ^ 



(13) 
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This shows that the NDC originates in the voltage- 
dependent occupation numbers and that a sufficiently 
strong occupation inversion (M12 < 0) and/or sufficiently 
fast reduction of the orbital polarisation (&M12/&V < 0) 
are necessary. We note that the second orbital is not oc- 
cupied for negative voltages, thus NDC is absent in this 
regime. For generic B and r NDC shows up in the vicin- 
ity of ni^ = n 2 -|- in the interval 5 < V < S + B. NDC 
is favoured for small magnetic fields, as the occupation 
inversion becomes less feasible with increasing Zeeman 
splitting. While the dot occupations do not exhibit re- 
markable renormalisation effects with respect to the PT 
results (see fig. @k), the conductance is strongly affected 
by both renormalisation effects and orbital PS contribu- 
tions. The NDC is strongly enhanced by the RG, and in 
contrast to the case of symmetric couplings (see fig. [3Ji) 
the orbital PS has opposite effects on the PT and RG 
results: while the orbital PS reduces the NDC in PT, 



it induces an increase in the RG results. Depending on 
the initial system parameters the NDC might even exclu- 
sively be generated during the RG flow. At low bias volt- 
age the effect of the orbital PS contributions is strongly 
reduced in PT, from the analytic calculations this effect 
can be shown to be proportional to (P ) 2 ~ 1/r 2 . In 
contrast, the corresponding RG results show a marked r- 
and ^-dependence of the orbital PS contributions. Con- 
cerning the cascade effect the onset of n^-f becomes more 
pronounced with stronger asymmetry. 

We stress that the mechanism leading to NDC in the 
present model is fundamentally different from the NDC 
observed in the interacting resonant level model [l8[ 
where it has its origin in a voltage-dependent renormali- 
sation of the hopping amplitudes towards smaller values. 

Finally we note that the temperature considered here 
is T = 0.01(5 > T K with the rough estimate T K ~ 
Dq exp[— 1/(NqJo)]. For higher temperatures T ~ 0.1 6 
thermal smearing inhibits a clear distinction of the dif- 
ferent transitions. The small absolute values of the con- 
ductance are due to the small values of the initial cou- 
plings NqJq ~ 0.01 required by our perturbative treat- 
ment. However, we expect the cascade effect as well 
as the appearance of NDC to be generic features of 
non-equilibrium transport through multi-orbital quan- 
tum dots. In particular, these features are expected 
to show up in the kind of transport experiments on 
molecules and nanowires as carried out in Refs. [2|- 
7] where the couplings to the leads will generically be 
strongly asymmetric and orbital dependent. 



CONCLUSION 



We studied a minimal model for a multi-orbital quan- 
tum dot exhibiting the cascade effect and NDC, namely 
a two-orbital Anderson model with asymmetric and 
orbital-dependent hopping amplitudes in a magnetic 
field. We considered the regime of strong Coulomb in- 
teractions where a SW transformation leads to an ef- 
fective Kondo model including orbital PS terms. The 
non-equilibrium transport through this system was stud- 
ied using a perturbative RG treatment with frequency- 
dependent couplings. In particular, the RG flow leads 
to a strong enhancement of the NDC. The observed 
strong influence of the orbital PS on the scaling of the 
Kondo couplings indicates an interaction-driven interplay 
of spin and orbital degrees of freedom even in absence of 
spin-orbit interactions or spin-dependent tunneling am- 
plitudes. The inclusion of additional spin-orbit inter- 
actions represents an interesting extension beyond the 
scope of the present analysis. 
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